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Abstract 

Heavy flavour decay to light hadrons is the key to understanding many aspects of the Standard 
Model from CP violation to strong dynamics. It is often presumed in line with the simple quark 
spectator model of D decay to Ktttt that the Ktt system has only / = 1/2. E791 have recently 
presented an analysis of their results on D + — > {K~tt + )tt + using a generalised isobar picture of two 
body interactions. While higher Ktt waves are described by sums of known resonances, the S'-wave 
amplitude and phase are determined bin-by-bin in Ktt mass. The phase variation is found not to be 
that of K~tt + elastic scattering. This hints at a different mixture of I = 1/2 and 1 = 3/2 S-wave 
interactions than in elastic scattering. Applying Watson's theorem to this generalised isobar model 
allows us to estimate the 1 = 3/2 Ktt S'-wave component. We indeed find that this is larger than 
in hadronic scattering or semileptonic processes. 
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Drawing quark line diagrams for heavy flavour decays provides a guide to the weak interaction 
dynamics that takes place. Thus as shown in Fig. la, in D decay to Ktttt, the c quark is 
seen to change into an s quark by emitting an off-shell W + that materialises as a tt + . The 
other quark in the D, for instance a d, acts merely as a spectator and with the s quark 
forms on the creation of a uu or dd a Ktt system. Since the isospin of the c and s quarks 
is zero, it is then natural to assume that the Ktt system has only the isospin of the u or d 
quark, namely J = 1/2. This presumption is in keeping with a simple isobar picture of this 
decay in which well-known 7 = 1/2 K* resonances, like the _fT*(890) or ^(1430), are readily 
seen. However, there are, of course, tttt resonances too, like the p, in other charge channels, 
that clearly indicate that the spectator picture is too simplistic. The W can be internally 
converted into hadrons as shown in Fig. lb, and then the Ktt system can have 1 = 3/2 
as well as / = 1/2. From wholly hadronic reactions we know the 1 = 3/2 interactions do 
not produce resonances, in keeping with the qq model of hadrons. Moreover, in low energy 
hadron scattering unitarity is a key constraint and restricts the magnitude of every partial 
wave amplitude to be bounded by unity. However, in the weak decay reaction we have no 
such constraint. Whilst in K decays we know there is the well-known AI = 1/2 rule, no such 
property is known for charmed decays. Consequently, we cannot presume that the 1=1/2 
component always dominate over 1 = 3/2 Ktt interactions, even though only the former is 
dominated by resonances. Here we estimate how big the fraction of / = 3/2 amplitude is in 
D decay. 



D + 




Figure 1: Quark line diagrams for .D-decays. (a) is the typical spectator model that underlies the 
idea that the resulting Ktt system has just I = 1/2. (b) involves internal TF-boson conversion. For 
semi-leptonic decays only the analogue of (a) is possible. This is shown as (c). 

This is made possible by a recent model-independent treatment of the E791 results on 
D — > Ktttt decay presented at La Thuile by Meadows pQ. There a Dalitz plot analysis is 
performed that isolates the phase and magnitude of the Ktt S-w&ve interaction bin by bin 
in Ktt mass. To explain how this allows us to estimate the 1 = 1/2 and 3/2 components, we 
begin with Ktt — > Ktt scattering. For this the S"-wave amplitude, in particular, is given by 

T(Ktt^Ktt;s) = d/aT^OO + C 3/2 T 3 /\s) , (1) 
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where s is the square of the Ktt mass, and the Cj are the appropriate Clebsch-Gordan 
coefficients depending on the charges of the Kn system. The hadronic amplitude T 1 with 
definite isospin has a magnitude and phase related by elastic unitarity for */s < (m^ + m^), 
the first strongly coupled threshold, and given by: 

T\s) = - sinS 1 exp(^ 7 ) , (2) 

where p = 2k/ y/s with k the Kir cm. 3-momentum. For the famous LASS experiment j2] 
the K~tt + — > K~n + amplitude is thus given by 



T(K-n + -> K-tt + ; s) = |- [sin<5 1/2 exp(^ 1/2 ) + \ sin 5 3/2 exp(^ 3/2 ) 



(3) 



These data, plus results on K~tt~ — ► K~n~ [3 , fix 5 1 for both / = 1/2, 3/2. The LASS 
group provide parametrisations 4 3 of these data, which are displayed as the dashed lines in 
the lower part of Fig. 2. 

Though these are valid in the experimental region LASS accessed, i.e. above 825 MeV, their 
extrapolations down to Kn threshold are known to be inconsistent with Chiral Perturbation 
Theory. Consequently, we consider a parametrisation [S] of the elastic phases, shown as the 
solid lines in Fig. 2, that is consistent with the calculations of Chiral Perturbation Theory jB] 
as shown in the dispersive analysis by Buttiker et al. [Zj These phases are the key ingredients 
in our discussion. 

We now turn to D decay, firstly in its semi-leptonic mode to Knfiu. The S'-wave amplitude 
in the Ktx channel is given by 

T{D - (Kn)w s) = F l J\s) + ^J\s) , (4) 

where the ^(s) are the semileptonic production amplitudes with I = 1/2,3/2. In this 
case the emitted W + is the source of the dilepton system, which of course, has no strong 
interaction with the Kn pair, as illustrated in Fig. lc. Now elastic unitarity requires that 
the phase of the Kir interaction with definite spin and isospin is the same as that in elastic 
Kn scattering, so that: 

F'M = I *S00 I exp [i5\s)] . (5) 

This is the famous final state interaction theorem of Watson |8,. Results from FOCUS [9 
confirm that this relationship holds and that the 1 = 3/2 component is small or negligible 
in this process. 

With this confirmation of Watson's theorem in heavy flavour decays, we consider the wholly 
hadronic channel: D + decay to K~7t + tt + . This has been analysed by E791 in a "new 
approach to the analysis of 3 body decays" pQ . The structure of the Dalitz plot has marked 
bands in both K~tt + mass combinations for the well known i^*(892) and iT|(1430) that 
feature so prominently in K~tt + — ► K~n + scattering |2j- This suggests the dominance 
of 2-body strong interactions. Consequently, E791 assume a generalised isobar picture in 
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Figure 2: In the upper and lower figures are the magnitudes and phases, respectively, of the 
1 = 1/2 and 3/2 amplitudes, T 1 , given by fits (black line) to the LASS results on K~ir + — * K~ir + 
scattering, shown as the data points. The fits shown are built from the 1 = 1/2 and 3/2 phases 5 1 
according to Eq. (3). The dashed lines are the LASS group's own parametrisation 0]. The solid 
lines [5] display a representation that not only describes the data, but extrapolates to threshold in 
keeping with one loop Chiral Perturbation Theory p]. 
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which genuine 3-body interactions are neglected. As is by now standard [TUJ E] , the decay 
matrix element is represented in the way summarised, for instance, by the CLEO-c group 
in Section II of Ref. PU| and so for the P-decay we consider here is then a sum of the 
amplitudes for the two possible combinations of K~tt + interactions with their appropriate 
spin and vertex factors. As in the usual isobar picture, all waves with J > are represented 
by a sum of resonance terms. This incorporates the known I = 1/2 resonances like the 
Pf*(892) K*(1A10) and K*(1680) for J = 1, for instance. With these resonances being well 
separated, their amplitude can be reasonably described by a sum of Breit-Wigner forms. 
Each partial wave has a production phase from the P7r — > K* coupling, which is assumed 
to be a constant as is also standard [HH [I]. Since the parameters of each resonance are 
those found in Kir scattering, e.g. by LASS [2], the phase variation of each Kit partial 
wave in P-decay is exactly that of Kit scattering in the elastic region, in agreement with 
the application of Watson's theorem to such an isobar picture [T2j. Since one overall phase 
is not determinable, the strong P-wave is taken as the reference wave with its production 
phase set to zero. 

What is new in the E791 analysis is the description of the Kir 5-wave. There the well- 
established wide Ifg (1430) appears, with perhaps an even broader k at lower mass. These are 
not simply describable by a sum of isolated Breit-Wigners, as the K~n + — > K~n + amplitude 
of Fig. 2 illustrates. Rather than enforce some ad hoc prescription of this key wave, E791 
represent this by a magnitude and phase in each bin of Kir mass in each K~tt + combination. 
The beauty of this analysis is that each 5-wave Kit mass band overlaps with a crossed Kir 
band in a P-wave and so the relative phase is determined, as well as the magnitude of the 
S'-wave. This provides as close to a model independent determination of the K~tt + S-wave 
interaction in this P-decay as is presently possible. The results are shown in Fig. 3. From 
the application of Watson's theorem to this generalised isobar picture |T2] , one would expect 
the phase variation of this S"-wave amplitude to follow that of Kit scattering in the region of 
elastic unitarity. As already mentioned (just before Eq. (2)), elastic unitarity is found [HI 12] 
to be very nearly exact upto Krf threshold, despite the fact that this is above the opening 
of the Kim and Kr\ channels, from which we infer these to be negligible below 1450 MeV. 
While the phase of the Kn S'-wave in P-decay and that of Kir scattering (compare the 
lower graphs of Figs. 2 and 3) both have an upward trend in the "elastic" region, they do 
not match to the precision expected. The reasons for this can be manifold: 

• the isobar assumption of only 2-body Kir interactions may not be true, 

• even if it is, then perhaps Watson's theorem does not apply and so the phase variation 
is not required to be the same, 

• even if Watson's theorem applies, should the phase variation of the K~7t + interactions 
be that of K~ir + — > K~ir + scattering or that of just its 1 = 1/2 component? 

These are all questions that have been raised at the BaBar Dalitz Workshop in December 
2004 ng. 
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Figure 3: The data are the magnitude and phase of the S-wave amplitude determined by the E791 
analysis of their .D-decay results. The "best" fit to these data with a quadratic form for the oti 
of Eq. (18) is the red curve labelled 2, also shown in Fig. 4. Curves 1 and 3 in blue and green 
respectively mark one standard deviation away from the optimal fit. 
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Since the generalised isobar picture of pairwise AB interactions defines a model for AB — > 
AB scattering, we adopt the view that this must be in keeping with data on AB — > AB 
scattering and the phase variation must agree in the region of elastic unitarity [12] . This has 
to be so, for the model to be consistent. Thus the Kir S-w&ve phase variation found by E791 
must agree with that for Kir — > Kit scattering. However, unitarity is only diagonalised by 
partial wave amplitudes with definite quantum numbers such as isospin, and not in general 
for individual charged channels, like K~n + . Consequently, it is the phase variation with 
definite isospin that should match. There is no reason that the combination of I = 1/2 and 
1 = 3/2 contributions determined by simple Clebsch-Gordan coefficients in K~tt + scattering 
(as in Eq. (3)) is that formed in D-decay, to which we already alluded in the introduction. 
The relationship between the phase variations then allows us to estimate the relative 1=1/2 
and 3/2 components of the E791 5-wave, as we shall shortly describe. Of course, there could 
also be I = 3/2 components in the higher partial waves too. These have been neglected in 
the E791 resonance-dominated description of these amplitudes. Nevertheless, it is known [3| 
that the J > 1 waves with I = 3/2 vary by less than 3° between Kn threshold and 1.8 GeV. 
Consequently, any such components should have only a tiny effect on our estimate. 

In keeping with the generalised isobar description that E791 adopt, the effect of the spectator 
pion is to produce an additional production phase, Pi, which we take to be a constant, as 
they do for all other waves. The K~tt + S-wave determined by E791 we call T , which is a 
function of s, the square of the Kit invariant mass, E. This is then given by 

Hs) = FlLlis) + = A(s) expMs)} , (6) 

where the amplitudes with definite quantum numbers are given by: 

^ I (s)had = I Fhad( S ) I eX P [ i§I ( S ) + iPl] > ( 7 ) 

in the region of elastic unitarity below 1450 MeV. The phases f3j reflect the structure of 
the complete set of quark line graphs of Figs, la, b in the strong coupling limit before final 
state interactions are included. Since different graphs contribute to each Kir isospin, the 
production phase @j depends on /. It is the magnitude A and phase, 0, of the total S'-wave 
Kit amplitude, displayed in Fig. 3, that E791 have determined pQ. Given this and the Ktc 
phases, S 1 shown in Fig. 2, the aim is to deduce the magnitude of the 1 = 1/2 and 3/2 Ktc 
amplitudes in D decay. 

Let us see how to do this. In vector terms it is like finding the vectors T 1 / 2 and JF 3 / 2 knowing 
only their sum A, i.e. 

A = F 1 ' 2 + T*l 2 . (8) 

Clearly there are an infinite number of vectors T Y that satisfies this. If the production phases, 
Pi, were zero (or otherwise determined), the solution would be simple, since we know from 
Kit elastic scattering what the phases 6i are at each energy in the elastic region, Fig. 2. 
Knowing the directions of the vectors, their magnitudes can easily be found. However, here 
we do not know in advance what the production phases are, but we do know that these 
are the same phases at every Kn mass, E = yfs. This sets the scene for determining the 
amplitudes. 
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Starting from 

A exp(#) = T 1 ' 2 + T 3 ' 2 , (9) 

where the T l are the complex I — 1/2, 3/2 amplitudes, it is straightforward to check that 
the solution is : 



tm(5^(E!)-^(E!)-p 1/2 + p 3/2 ) " I >\ 

sin f<jV2(^) + -0(E)) r , , 

sm (d 1 /^) - <Jd/2(£!) + /3 1/2 - /3 3/2 J 



The phase difference, 5^ 2 (E) - 5 3 / 2 (E), varies from to ~ 180° in the elastic region below 
Kv( threshold, and so we see that the denominator in each expression, namely 

sin (S^E) - 5 3 ' 2 {E) + p 1/2 - p 3/2 ) (12) 

will inevitably vanish at at least one energy in this region. Let this reference energy be E r . 
Since the amplitudes T l are finite at real energies, the numerator must vanish at this same 
reference energy. This fixes the production angles. Thus 

A/2 = 4>{E r ) - 5 1/2 (E r ) + mn (13) 
A/2 = 4>{E r ) ~ S 3/2 (E r ) + TVK , (14) 

where m, n are integers (including zero). Consequently, we can write the individual isospin 
amplitudes in terms of just one unknown parameter, the reference energy E r , as: 

= A sin(m-m)-^ 2 {E) + 5 3 l 2 {E r )) 



sin (5 1 /2( jB ) _ 5 3/2( E ) _ sV*(E r ) + 5^ 2 {E r )) 
x exp [i{5 1/2 {E) - 5 1/2 (E r ) + <f>(E r ))] (15) 

sin (^(E)-^/ 2 (E r )- <P(E) + 0(E T .)) 

M Sin (5V2( jE ;) _ §3/2(E) - 8^{E r ) + £3/2(£ r )) 

x exp ^(5 3 / 2 (E) - 5 3/2 (E r ) + 0(£ r ))] . (16) 

For each value of E r we can then determine the amplitudes. Consequently, we have a 
continuous range of possible amplitudes. Within this set, there are only a small group 
that are physically (as opposed to mathematically) allowed. Watson's theorem is ensured 
by the phase variation of each of JF 1 / 2 and T 3 ! 2 being given by the elastic phases of the 
corresponding elastic scattering amplitudes T 1 . However, the dynamics of the final state 
interactions requires that the amplitudes T 1 should be smoothly connected to the scattering 
amplitudes T 1 through coupling functions aj(E). A consequence of this is that any resonance 
poles in Kn scattering appear in the decay process with the same mass and width. The 
functions SJ fix the coupling to the decay channel. Thus 

F\E) = aj(E) T\E) exp(i/3j) . (17) 
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Figure 4: The moduli of the 1 = 1/2 and 3/2 Kir components of the -D-decay amplitude determined 
from the E791 analysis using Eqs. (15,16) with E r ~ 800 MeV. The error bars shown are fixed by 
the uncertainties given by E791. They are inevitably largest in the region of E ~ E r . The fit 
given by representing the coupling functions ckj as quadratics in s is shown. 
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The dynamics of the pseudoscalar interactions imposes Adler zeros in each of the Ktt scat- 
tering amplitudes with / = 1/2, 3/2 at s = Sj . Such zeros may appear in the decay 
process but not necessarily at the same position as in elastic scattering. We take account 
of this by defining reduced elastic scattering amplitudes with the Adler zero divided out. 
We then specify new coupling functions, oti which have to be smooth since they contain no 
explicit s— channel dynamics. Thus we have 

F\E) = ai (E) exp(^) , (18) 

S S j 

where the coupling functions aj are representable by low order polynomials in s = E 2 . 

We therefore determine the coupling functions ai(E) for each choice of E r and select those 
that are representable by low order polynomials, while still providing an accurate fit to the 
amplitude determined by the E791 analysis. With 26 values of magnitude and 26 phases in 
the elastic region below Krj' threshold, this is a severe constraint. 

The "best" fit is found with E r ~ 800 MeV, for which the production phases are @i/ 2 = 72° 
and Ps/2 = —73°. The moduli of the / = 1/2 and 3/2 amplitudes, using Eqs. (15,16), are 
shown in Fig. 4. As can be seen these are accurately determined away from E r , where the 
errors are very small. These show that the ckj of Eq. (18) are well represented by quadratic 
polynomials in the square of the Kn mass, s. The precise results above 1 GeV fix the 
continuation to lower Kn mass as illustrated in Fig. 3. The requirement that the Z)-decay 
amplitude of Eq. (6) should be well described is illustrated by the (red) lines, labelled 2 in 
Fig. 3. Allowing a one standard deviation change by varying E r alters the 1 = 1/2 and 
3/2 amplitudes, while still fitting the E791 magnitudes and phase. This result is shown in 
Figs. 3, 5. The lines 1, 3 (in blue and green, respectively) in Fig. 3 show the corresponding 
small difference in the description of the E791 amplitude, while Fig. 5 shows the range of 
variation of the 1=1/2 and 3/2 then permitted. The bands reflect a change of E r from 750 
to 950 MeV, or equivalently a range of 

1/2 = 65° -> 85° , /3 3/2 = -86° -> -38° , (19) 

which are strongly correlated. 

We, of course, only know that the coupling functions of Eq. (18) should be smooth 
functions of Ktt mass since they do not contain direct Ktt dynamics. A constant or linear 
function of s does not give an acceptable fit in terms of \ 2 for any value of E r . A quadratic is 
the lowest order polynomial to give fits of acceptable confidence, as seen from the curvature 
of the magnitudes in Fig. 4 determined by the data-points with small error bars above 1.1 
GeV, as shown in Fig. 3. Adding higher order terms in the polynomial representation of the 
ai does not significantly improve the confidence level. 

We see by comparing the upper plot of Fig. 2 with Fig. 5 that the 1=1/2 and 3/2 Ktt 
components of D decay are quite different from those of elastic scattering. The 1 = 3/2 
component is more than 50% of the 1 = 1/2 above 1.1 GeV. This is surprising because the 
/ = 1/2 amplitude contains the broad _K"q(1430), while the I = 3/2 is entirely non-resonant. 



10 



12 



I 



10 - 



2 - 



4 - 



8 - 



6 - 








0.6 



0.7 



0.8 



0.9 



1.1 



1.2 



1.3 



1.4 



M Kn (GeV) 



Figure 5: The red (solid) lines give the magnitudes of the I = 1/2 and 3/2 components of the 
D-decay amplitude determined by the "best" fit to the data of Fig. 3. The bands delineate the 
ranges set by one standard deviation variation from this "best" fit. 

Consequently, any analysis of the D-decay data which neglects the 7 = 3/2 amplitude must 
place this component incorrectly in other contributions and so lead to false conclusions about 
resonant branching fractions. Neither isospin component appears to contain an Adler zero 
and hence they grow towards threshold. At low mass we see the I — 1/2 amplitude is 
dominant, whether this is because there is a near threshold k resonance requires a far better 
determination of the amplitude and phases. 

An important ingredient in fixing the isospin components of the D-decay amplitude are 
the phases of Kir scattering in the elastic region. Here these have been set by the LASS 
experiment above 825 MeV and their continuation down to threshold by the predictions of 
Chiral Perturbation Theory. However, increased statistics on D — > Ktttt decay would reduce 
the error bars on the amplitude and phases in Fig. 3. These might well become sufficiently 
small that an acceptable fit cannot be obtained without changing the K-k phases 5 1 away 
from their LASS results. This would mean that the D-decay results could indeed increase 
the precision within which the Kit elastic phases are known. Such improvement in D-decay 
statistics would be particularly welcome below 1 GeV down to threshold, where the current 
errors are sizeable — Figs. 4, 3. There it holds out the prospect of revealing whether there 
is indeed a low mass k or not. Such a state can only be exposed by analytic continuation of 
the appropriate amplitudes into the complex s-plane to deduce whether a pole exists or not. 
This requires far greater precision than the E791 data and its analysis we have considered 
here. D-decay results from D-factories may make this possible. The present analysis points 
the way. 
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